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The relativistic theory for the electric dipole moment (EDM) of paramagnetic atoms arising 
from the electric dipole moment of the electron is presented. A novel approach using the 
relativistic coupled-cluster method that incorporates the residual Coulomb interaction to 
all orders and a weak parity and time-reversal violating interaction to one order has been 
employed in Fr to obtain the enhancement of the EDM of that atom compared to the EDM 
of the electron. Trends of the different correlation effects and leading contributions from 
different physical states are discussed. Our result in combination with that of the Fr EDM 
that is currently in progress, has the potential to probe the validity of the Standard Model 
(SM) of elementary particle physics. 

PACS numbers: 

I. INTRODUCTION 

The Standard Model (SM) of elementary particle physics given by Glashow, Weinberg and Salam 
described by the SU{3) x SU{2) x U{1) gauge group satisfactorily explains most of the phenomena 
observed so far on the fundamental interactions [l|, [3, Q]. However, this model is widely believed to 
be incomplete. Indeed a question of great importance is whether there is any new physics beyond 
the highly celebrated SM of elementary particles. Any kind of experimental evidence supporting 
the extensions of the SM like the left-right symmetric model 4|, the multi-Higgs model Q] or 
supersymmetry models 6!] will have a significant impact on our understanding of the fundamental 
forces between the elementary particles in nature. 
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Detailed studies of the fundamental interactions have a bearing on chemical and biological 
processes as well. For example, the violation of parity is believed to be one of the fundamental 
reasons behind the "homochirality" of bio- molecules [t]- In the same vein, any time reversal 
symmetry (T) violation will lead to a certain asymmetry in chemical phenomena whose importance 
will depend on the actual magnitude of the asymmetry observed in the experiments. A case in 
point is the emergence of enantiomorphic excess where none is expected from the electromagnetic 
interaction Q]. One can in fact, distinguish between the chirality brought out by T-invariant and T- 
noninvariant enantiomorphism, introducing the concept of "false chirality" in the latter situation. 
Clearly if the predictions by the SM and those beyond, relating to certain symmetry violations 
are observed, they will have profound chemical and biological implications. Additionally, if the 
predictions from the SM and some of its extensions differ, then this will lead to different expected 
magnitudes of the effect. This may bring observational efficiency to the limit achievable currently. 
Fixing the limit of observability then depends on the accuracy of electronic structure theories [8]. 
This is why studies on the predictions as defined by the SM and beyond are of immediate interest 
to chemists. 

The EDM (D) associated with the nondegenerate state of a physical system like an atom can 
be shown to be proportional to the angular momentum (J) of that on the basis of the Wigner- 



Eckart theorem 



^. It can be proved that the existence of a nonzero electric dipole moment (EDM) 



which is the expectation value of operator D in a 



given atomic state would imply the violations 



of the T invariance and parity (V) symmetries [lO|, The SM which conserves the combined 
transformations of charge conjugation (C), V and T symmetries (known as the CPT theorem jl^). 
predicts V and CV violations in weak interactions. A possible observation of a nonzero EDM in 
an atom would be a direct signature of T violation. An EDM has not been observed so far for 
whether in elementary particles or composite systems. The SM predicts an upper limit on the 



electron EDM of the order of 10 ^^e ■ 



cm 
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l4\ . while other suggested models predict this value 



almost ten orders of magnitude larger. Therefore, measurements of EDM in atomic systems are 
necessary to test the SM. Thus, an unambiguous observation of a nonzero EDM of any fundamental 
particle including the electron will undoubtedly unveil a new arena of physics beyond the SM. A 
direct measurement of the electron EDM is not possible using accelerator approaches due to the 
charge of the electron, however it could be obtained by combining the results measurements and 
calculations of atomic EDMs using the table-top experiments. 

In 1963, Schiff [15] argued that the permanent EDM of an atom vanishes under the assumption 
that the constituents of non-relativistic point-like systems interact only electrostatically. It was. 



however, put forth by Salpeter [3] and Sachs and Scwebel [l^] as early as 1958-1959 that when 
relativistic effects are included it is possible for an EDM on the electron to give rise to an EDM 



111 ] that by introducing 



on the atom as a whole. Later, it was rigorously proved by Sandars 
the correct relativistic form of the interaction Hamiltonian into the Dirac equation in a Lorentz 
covariant manner, as emphasized by Salpeter, one will get non-zero EDM for the atom when an 
electron is assumed to have intrinsic EDM. 

In principle, an atom being a composite many-body system will have many sources of intrinsic 
EDMs arising from its constituent particles and their V, T violating electron- nucleus interactions. 
At the elementary particle level; electrons may possess intrinsic EDMs, describing by a coupling 
constant de which may directly contribute to the atomic EDM. The electron EDM can also interact 
with the central electrostatic field of an atom and can produce an atomic EDM. The interactions 
between electrons - quarks may manifest at different levels as electrons - nucleons to electrons - 
nucleus interactions. The V and T violating electron- nucleus interactions can either be scalar- 
pseudo-scalar (S-PS), with a coupling constant Cg or a tensor-pseudo-tensor (T-PT) interaction 
describing with a coupling constant Ct- 

The paramagnetic (open-shell) atoms will in general be sensitive to the contribution from the 
intrinsic EDM of the electrons and their and T-violating S-PS interactions with the nucleus 
whereas, the EDM in diamagnetic (closed-shell) atoms will arise from the EDM of the nucleus 
and their "P, and T-violating T-PT interactions with the electrons. This can be understood as the 
electron EDM is a spin dependent property, hence it contributes to open-shell atoms because of 
the unpaired valence electron. However, for a closed-shell atom because of the Fermi exclusion 
principle, it will add up to zero if hyperfine interactions are excluded. These V- and T-violating 
electron-nucleus interactions will provide a useful tool in understanding the CP violation from 
the semi-leptonic sectors. Using the knowledge of CP- violation obtained in these atomic EDM 
experiments one can also constrain different models of CP-violation. 

It was shown by Sandars that by choosing suitable atoms and favorable electronic states one can 
get an enhanced EDM for an atom which may even be a few orders of magnitude larger than that 



of the free electron [13] • In particular, he demonstrated by carrying out relativistic calculations 
that the atomic EDMs of thallium (Tl) and caesium (Cs) are two orders of magnitude larger than 
the EDM of the electron. It was this important result which provided the impetus for the first 
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2l[ |. It was 



generation of atomic EDM experiments that were carried out in the 1960s 
realized quickly that the EDM enhancement factor R defined as the ratio of atomic EDM to the 
electron EDM increases with increase in the nuclear charge (Z) and also if there are close lying 
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states of opposite parity since it is inversely proportional to the difference in their near degenerate 
energy levels. Hence, heavy rare-earth atoms with anomalously close energy levels of opposite 
parity, have large EDM enhancement factors. These factors are also fairly large for heavy alkali 
atoms. The enhancement factor is proportional to 22 



^ " J(J+1/2)(J + 1) ' 
where a is the fine structure constant. This formula for order of magnitude estimate of EDM 
enhancement factors illustrates the dependence on nuclear charge Z and the angular momentum 
J which implies that R is large for high Z and low J. 

The ground state EDMs of heavy neutral alkali atoms are of considerable interest to experi- 
mentalists, because: (i) the EDMs in these cases are several orders of magnitude larger than that 
of the electron, unlike in the hydrogen atom where the enhancement factor is large only for the 
2s excited state, however, it can be easily perturbed by external electric fields and is therefore 
not amenable to sensitive experiments, (ii) the availability of commercial lasers whose operating 
frequencies match with those of the resonant energies needed for causing transitions between the 
low lying levels and (iii) the large polarizabilities of these atoms 

The EDM of an atom or any other neutral particle is determined experimentally by applying 
an external static electric field to the atom and measuring its shift in energy that results from the 
interaction of the EDM with the electric field. Consider an atom which has a permanent EDM as 
well as a magnetic dipole moment. In the presence of a static electric field E and a magnetic field 

the interaction Hamiltonian is given by 

Rint = -D.E - p.B, (2) 

where D and /2 are, respectively, the electric and magnetic dipole moment operators. 

The application of the external fields leads to a precession of the atom. The precession (Larmor) 
frequency is primarily due to the magnetic dipole moment, but there is also a small contribution 
from the EDM. The observable in an EDM experiment is the difference in the Larmor frequencies 
corresponding to parallel and antiparallel configurations of E and B - reversal of E relative to B. 
This change in frequency is 

'^DE 

^^E = (3) 

It corresponds to Auje = W~^Hz for D~ 10~^^e — cm and E = WKV/cm. This frequency shift 
corresponds to a magnetic field of lO^^G for a diamagnetic atom and 10^^^ G for a paramagnetic 
atom. 
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One of the most important systematic errors in the EDM experiment is the magnetic field that 
is produced by the motion of the atoms. This field to first order in ^ is given in the moving frame 
of the atoms as 

B^ = -xE. (4) 

c 

For V = 300 m/s and E = 10 KV/cm, the motional magnetic field Bm = SxlO^^G. This magnetic 
field can give rise to a frequency shift which can mimic an EDM. In the Tl experiment, two 
counter-propagating beams are used to minimize this effect 



2j] . EDM experiments using optically 



pumped atoms in a cell have a zero average velocity and are therefore not affected very much by 
the motional magnetic field 



25( 1 . Both the beam and the cell experiments have their advantages 
and disadvantages. While it is possible to apply larger electric fields in the beam experiments, the 
coherence times are longer in the cell experiments. The motional magnetic fields often limit the 
sensitivity of the former, while leakage currents give rise to systematic errors in the latter and they 
cannot be estimated easily. 

EDM experiments based on laser cooled and trapped atoms in principle have the advantages of 



both the beam and the cell experiments 



26 



27 



28|, 



30l |. In these experiments, one can apply 



large electric fields and the coherence times are long. The leakage current problem can be overcome 
using a suitable configuration for the laser trap. The systematic error due to the motional magnetic 
field is virtually non-existent because of the extremely low average velocity of the cold atoms. 

The procedure of the atomic EDM measurements with cold atoms is as follows: First, a fast 
atomic beam from a hot oven is slowed down by using the Zeeman-tuning method. After this 
pre-cooling stage, the atoms are trapped and cooled by a magneto-optical trap (MOT). A high 
density and large number of atoms is then loaded into the MOT within several seconds. Then the 
atomic beam and the magnetic field for the MOT is switched off, the detuning of the trapping laser 
is increased and its intensity reduced. This results in further cooling to the micro Kelvin region 
by the polarization-gradient method. The next step is to perform optical pumping to polarize 
the nuclear spin by the application of a circularly polarized resonant light pulse after completely 
switching off the laser fields. Finally, a high power laser for a far blue-detuned dipole forced trap 
and also a high static electric field for the EDM measurement is switched on. A probe laser beam 
will be used to measure the Larmor precession frequency. The loading and measurement procedure 
is repeated many times to reduce the statistical uncertainty. 

Atomic theory is needed in combination with experiments to extract a variety of T violating 
coupling constants. Following a series of insightful calculations of atomic EDMs by Sandars in the 
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1960s and 1970s based on the relativistic central field potential ll|, a number of relativistic many- 
)ody calculations on atoms of experimental interest have been carried out in the last two decades 
3ll |. Ingenious ex per iments were initiated in the 1980s to observe the EDMs of mercury (Hg), 



Cs and Tl 



25 



32| . Subsequently systematic errors have been improved in these experiments 
and today the most accurate results are available in Hg for diamagnetic atoms and Tl for 
paramagnetic atoms [34]. Polar molecules seem to be better candidates than heavy atoms for 



observing EDMs arising from the electron EDM and there have 
the observable in some of the molecular EDM experiments 



)een recent attempts to calculate 



35|, ISg, 137[ . However, these calculations 



are in their infancy. In contrast, it is possible to perform very accurate calculations on alkali atoms. 
Therefore, to shed light on physics beyond the SM or in fact to observe a nonzero atomic EDM, the 
heaviest of the alkali atoms, francium (Fr) seems to be another suitable candidate for its relatively 
simple spectroscopic levels and and large enhancement factor. A number of spectroscopic analysis 
on Fr has been carried out at SUNY Stony Brook using laser trapping technique 
been proposed both for EDM 



and it has 



39( 1 and atomic parity violation (APV) [40] measurements. In fact, 



there is a preliminary calculation available for the enhancement factor of Fr using radial integrals 



for a few matrix elements 



411 ]. In the present work, we will demonstrate an approach using the 



relativistic coupled-cluster method that includes Coulomb interaction among electrons to all orders 
and V and T violating electron-nucleus pseudoscalar interaction due to electron EDM up to first 
order but involving contributions from continuum and all possible single and double excited states 
to calculate R in Fr. In this procedure, we will also demonstrate contribution from important 
intermediate states and highlight the role of various electron correlation effects. 

II. THEORETICAL STUDIES OF ATOMIC EDM 
A. General features 



The three discrete symmetries - charge conjugation (C), space inversion or parity (V) and 



121 ]. According to this theorem. 



time-reversal (T) are related by the well-known CPT theorem 
while a physical system described by a local field theory may violate any one of these symmetries 
independently, it is invariant under the combined operation of all three of them. 

The combined operation of C and V (CV) was thought to be a good symmetry till 1964 when 



Christenson et al observed its violation in the decay of Kq meson 



42| . Using the CPT theorem, one 



arrives at the conclusion that this observation implies T violation. An observation of an electric 
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J 




FIG. 1: Violation of T and V symmetries is necessary for a non-zero electric dipole moment. 

dipole moment (EDM) of a non-degenerate physical system would be a direct signature of the 
violations of T as well as V symmetries [43]. Consider a physical system with non-zero angular 
momentum that has a permanent EDM, D. The vector D like the magnetic dipole moment will 
be parallel or antiparallel to J, the angular momentum of the system [18]. Under the action of T, 
J — J and D ^ D. Since D is proportional to J, T invariance implies D = 0. If a transformation 
is carried out under V, D ^ —D but J ^ J, and this leads to D = for V invariance. One therefore 
concludes that T and V must be violated independently for the existence of a non-zero EDM. The 
above arguments are illustrated in Fig. [TJ 

B. Atomic EDM from the electron EDM 

In a manner analogous to the anomalous magnetic moment, one can introduce an electric dipole 
moment for a single electron in an electromagnetic field in a Lorentz covariant manner into the 
Dirac equation. The only difference is that it contains a pseudo-scalar Dirac operator 75. 
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The pseudo-scalar perturbation Hamiltonian for the intrinsic EDM of an electron then reads, 
Hedm = (^757^7;.^) F^^ = -^ {^l^^r^u'f) F^^ (5) 

where dg is the intrinsic EDM of the electron. 

In the Pauli approximation, Eq. ([5]) reduces, in atomic units, to 

Hedm = -deP [a-E + ia-H^ (6) 

where a, (3 and a are the Dirac matrices, E and H are the total electric and magnetic fields at the 
site of the electron. The second term in the rhs of Eq. ^ is much weaker than the first term and 
hence, we will consider only the latter in the calculation of EDMs of paramagnetic atoms. The 
non-relativistic limit of this interaction is — deff ■ E. This form of the interaction produces a zero 
atomic EDM thus as we show below the EDM of a paramagnetic atom due to the intrinsic 
EDM of the electron is entirely a relativistic effect. Therefore, it is necessary to use a relativistic 
many-body theory to determine its size. The atomic EDM arising from the intrinsic EDM of the 
electron is discussed in detail below. 



C. Method of calculations 



The total Hamiltonian for a many-electron atom, in the absence of any external field, when the 
electron possesses an intrinsic EDM is given by, 

H = Hq + Hedm (7) 

where, Hq is the atomic Hamiltonian given by 

Ho = Y, [ca -Pi + iP- l)c' + Vnuciri)] + — (8) 
t i>j ■' 

due to the electromagnetic interaction and Hedm = —deYli P^i ' -^i"* with the internal electric 

field, = -V Vnuciri) + Y.i>j y- > exerted by the nucleus {Vnuc{ri)) and other electrons (^). 

The application of external electric field also induces EDM to the atom. Hence, the total 

perturbed Hamiltonian H^^^ in the presence of an external electric field is given by, 

H^^^ = Hedm- deY,P^i-E - eY,n-E (9) 

i i 

In the time-independent perturbation theory, the first order shift in energy is given by, 

^(1) ^ /^(0)|^(l) |^(0)\ 

= _ (.I/W|4/3a,.^|^W)-(^W|er-;.^|^W)] (10) 
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Assuming that the apphed field is in the positive z-direction, 

3 

Noting the fact that, the operators in the first and third terms are odd under parity and the 
expectation value of odd parity operator vanishes in a state of definite parity, the only non- vanishing 
term (corresponds to the even parity operator) is given by, 

= -d,{^t^\(3Y:^a,,,\^t^)\E\ (12) 
i 

As the strength of the perturbation is sufficiently weak, we consider only up to the first-order 
perturbation in the wavefunction. The first-order perturbed wavefunction \^^) is given by. 



The second order shift in energy due to EDM as a perturbation is given by. 



™ ' 2^ p(0) _ p(0) ^^''> 



p(o) _ p(o) ^ ' 

On inserting the expression for H^^^ in the above equation and discarding the terms containing 
de^ and E"^ while expanding, we get those terms which are first order in perturbation and linearly 
proportional to the applied uniform electric field as below; 

f(o) _ f(o) 

+ ^(o)_^(o) jl^l (15) 

Thus, the total shift in energy would be. 



Ei^^ = { E 



f{0) p{0) 



+ E ^(o)_^(o) (16) 
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The coefficient of tlie linear sliift in energy due to the apphed electric field is the total electric 
dipole moment of an atom i.e, e = —D ■ E. Thus, the total EDM of an atom is given by, 



{D) = |4(^^^)|/?E^^.ii^-) - E 



p(0) p(0) 

2-^ p(0) _ p{0) f ^^') 

The total atomic EDM given by Eq. ()17p can be further simplified to the effective one-electron 
form given by, 



3 



O- ^ r /vI/(°)|/^'v5r.2|VI/(°)\/vI/W|? 1x1/(0) 

Zicde I v-^ \ Pjl^n ){^n \ Zj \ ^ m 



E 



p(0) p(0) 



+ h.c. (18) 



Here, /? and 7^ are the Dirac matrices, p is the 3-momentum, h (= ^) is the modified Planck's 
constant, E^^^ and are the zeroth order energies of the states m and n respectively, and the 



abbreviation h.c. stands for hermitian conjugate. Eq. (llSp shows that Hedm in Eq. (llTh reduces 
to an effective one-body term given by H'^^j^ = 2icde f^J^p]- 

^iIdm responsible for mixing atomic states of opposite parities. Its strength is sufficiently 

weak for it to be considered as a first-order perturbation. It is, therefore, possible to write the 
m'th state atomic wavefunction as 

|^^) = |*W)+4|*«). (19) 

Therefore, we have 

I / 

where D is the electric dipole (El) operator. 

The above first order perturbed wave functions due to H'^Jl/j^i can be calculated by summing 
over a few important intermediate states. However, the accuracy of this approach is rather limited. 
We have developed an approach based on the relativistic coupled-cluster (RCC) theory that can 
overcome the sum-over-states approach by directly solving the first order perturbed equation 

{H^'^ - E^^^)\^g^) = (41^ - H^JL)\'^t^), (21) 

where E^i^ vanishes since H'^j^j^,^ is an odd parity operator. We present below the formulation of 
this problem based on RCC theory. 
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D. RCC theory of atomic EDM 

Using RCC theory, the atomic wavefunction |^'m^) for a single valence (m) open-shell system is 



given by 



\^^^^)=e^'"'{l + S^^^}\^m), (22) 



where we define \^rn) = aJnl'^o)) with |$o) as the Dirac-Fock (DF) state for the closed-shell system. 
The curly bracket in the above expression represents normal order form. 

In the single and double excitations approximation coupled-cluster (CCSD) method, we have 

c{0) _ o(0) , r,{0) (^o) 

where T^^^ and Tg'^'' are the single and double particle-hole excitation operators for the core elec- 
trons and and are the single and double excitation operators for the valence electron, 
respectively. The amplitudes corresponding to these operators can be determined by solving the 
relativistic coupled-cluster singles and doubles equations. 

These amplitudes are obtained by solving the following equations 

(^.^|{^}|$0> = 6o,LAE,orr (24) 
= -{^i\{He^}\^m) - {^i\SL,rr. + S„|$™)IP, (25) 

where the superscript L(= 1,2) represents excited states from the corresponding DF states, the 
wide hat symbol denotes connected terms, AEcor and IP represents the correlation energy and 
ionization potential of the valence electron 'm', respectively. We consider the Dirac-Coulomb (DC) 
Hamiltonian in our calculation given by 

H = [Ho] + [v;,] 

N 

= K [c« ■P^ + {P- l)c' + Vnucir,) + C/(r,)] A+ 

i 

N N 

+E ^^-^^ - E ^tUinW], (26) 

i>j ^-^ i 

where Hq is the DF Hamiltonian, V^g is the Coulomb residual term in atomic units that is neglected 
in the DF calculation and A"*" are the projection operators on to the positive-energy states of the 
Dirac Hamiltonian in the nuclear (Vnudri)) and DF {U{ri)) potentials. 
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The most important triple excitations have been considered by constructing triple excitation 
operators [46|] 



1/ T-(0) , c(0) 
e?) + — Eg — 

where e's are the orbital energies. The above operators are used to construct single and double 
open-shell cluster amplitudes by connecting further with the CCSD operators and evaluating con- 
tributions to IP which are also involved in the amplitude determining equations. Finally, these 
amplitudes are solved self-consistently. 

In the presence of the EDM, the exact atomic wavefunction can be written as 

=e^{l + S^}|$„), (28) 

where the cluster amplitudes are given by 

Sm — S^yi^ -\- d^s!^\ (29) 

r(i) and are the first order in de corrections to the cluster operators T^'^^ and Sm^ , respectively. 
In the CCSD method, we have 

s^r^^ = + ^S- (30) 

The amplitudes of these operators are solved, keeping up to linear in de, by the following 
equations 



($^|<)ra)+i/^^^,|$o) = o, (31) 



{^tlHPs^r^^ + i/ff)r(i) + i/S)r(i)5L°) + H'Jl,, + H'Jl,,S^^^\'^^) = -($^|S«|$^)IP(32) 

where Hamiltonian operators with overline as defined as Hj\f = e~^^°\H}e^'°\ These are com- 
puted after determining T^^\ 

After solving the CCSD amplitudes we determine (D) by 

„ (D) <cI>^|{g^^^V(l + rW^)5^°^^}e^'°''l)e^'°'{(l + rW)gi°^+5£^ 

de (1 + A^r) 

In the above expression, we define D^^^ = e^^°^ De^^'^^ and Nm^ = {^rn.\Sm^ e'^'°' 6^^°'' Sm \^m) for 



the valence electron 'm'. These terms are evaluated using the generalized Wick's theorem 
constructing effective one-body, two-body etc. terms. 



m by 
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As explained in our earlier work on EDM arising from scalar-pseduoscalar interaction in Cs and 
Tl 0, the important RCC terms in the above expression are Z^WT^^^^ ^^^[11 and ^DWsJ^- The 
first term corresponds to the core-correlation effects whereas the second and third terms correspond 
to the valence correlation effects. Core-polarization and pair-correlation effects arising through the 
singly excited states are considered through and important core-polarization effects from 

the doubly excited states are taken into account through D^^^S^^. Important excited states that 
contribute significantly through these terms are given below. 



III. CONSTRUCTION OF BASIS FUNCTIONS 



The EDM interactions for both closed and open shell atoms are sensitive to the nuclear region. 
Therefore, the Gaussian type orbitals (GTOs) that produce good wave functions in the nuclear 
region can be used to calculate accurate {D). For atomic wavefunction calculations, they are given 



by 



i^.,fc(r) = r'=e-"''-', (34) 
where A; = 0, 1, .. for s,p,.. type orbital symmetries, respectively. For the exponents, we have used 

ai = ao/3'~^. (35) 
We consider excitations from all occupied orbitals (holes)in both the DF and RCC calculations. 



491 ]. The finite size of the 



The orbitals are generated on a grid as in the numerical code GRASP 
nucleus has been accounted for by considering a two-parameter Fermi nuclear charge distribution 
approximation given by 

P = 7^^- (36) 

^ 1 + e(r-c)/a ^ ^ 



We use values of po) c and a as given by Parpia and Mohanty 



IV. RESULTS AND DISCUSSIONS 



In Table [U we present our enhancement factor calculations to the ground state of Fr. We obtain 
the final result R = 894.93 which agrees fairly well with Byrnes et al. who have obtained 910(46) 



(error bar is quoted as ~ 5%) [4]|. However, this is just a coincidence. Byrnes et al. have used a 
sum-over-states approach and have only considered singly excited valence states as the intermediate 
states. They have therefore not taken into account pure core correlation and contributions from 
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TABLE I: Contributions from important RCC terms to the R — Da/d^ calculations of the ground state in Fr. 
c.c represents conjugate terms. Here contributions given by Norm and Others correspond to normalization 
corrections higher order terms those are not mentioned in this table, respectively. 



RCC terms 




R = Da/de 


Dirac-Fock (core) 


25.77 


Dirac-Fock (virtual) 


695.44 


Dt'i^^ + c.c. 




43.39 






1000.19 






-64.94 




- c.c. 

- c.c. 


-18.07 
-59.18 




- c.c. 

- c.c. 


-2.80 
19.26 


Norm 




-24.42 


Others 




1.51 


Total 




894.93 



the doubly excited states. In their calculation, the 7si/2 ~ ^Pij2 El matrix element was taken 
from experiment and only 7pi/2 ~ lOPi/2 discrete states were used in the calculation and the 
contributions of the continuum states were included in an approximate way. Another limitation in 
the calculation of Byrnes et al is in their treatment of the internal electric field in their parity and 
time-reversal violating Hamiltonian. They have considered only the potential that an electron sees 
from the nucleus but not the other electrons of the atom. We have included contributions from 
all the core-electrons and all possible single and double excited states through the RCC method in 
our calculation. 

We have also explicitly given the DF contributions from core and virtual orbitals in Table [H It 
is clear from this that the contributions due to core orbitals are not small in the present system. 
Our DT^^ + c.c. RCC terms at the lowest order correspond to the above DF core contributions. 
Comparing contributions at the DF and the RCC levels, it shows that the all order core correlation 
effects are almost twice as much as the DF contribution. Similarly, our + c.c. RCC terms 

at the lowest order corresponds to the DF results only due to the virtual orbitals. From Table [U 



it is also clear that the all order correlation effects are significant. 

TABLE II: Break down of DT^^"^ +c.c. contributions using DF reduced El matrix elements from core orbitals. 







{7si/2\\D\\npi/2)DF 


{npi/2\Ti^^\7si/2) 


R 


n 


= 2 


-0.15 X 10-3 


1.55 


0.19 X 10-3 


n 


= 3 


0.88 X 10^3 


-3.35 


0.24 X 10-2 


n 


= 4 


-0.52 X 10-2 


7.56 


0.32 X 10-1 


n 


= 5 


0.38 X 10-1 


-20.31 


0.63 


n 


= 6 


0.54 


-100.01 


44.16 



Contributions from the doubly excited states arise through the operator. It can be seen 
from the above table that these states contribute to about 7% iof the final result, but with opposite 
sign. Therefore, the error bar quoted (which is ~ 5%) by Byrnes et al. without considering the 
doubly excited states contributions does not seem to be correct. Again, the normalization correction 
(Norm) in our calculation is also about 3%, which is a significant contribution, and this appears 
to be missing in the work of Byrnes et al. 

The trends exhibited by the correlation effects in the present work are similar to those in the 



case of rubidium (Rb) 



51[, but the amount of core correlation effects are substantially larger for 



Fr. Interestingly, we observe that the correlation contributions for Rb [51], Cs [51] and Fr are 
about 24%, 22% and 20%, respectively. This is because of the fact that as the size of the system 
increases contributions from the doubly excited states increases with opposite sign and hence there 
are strong cancellations in the heavy systems. 

In Table [Til present the core orbital contributions combining the reduced matrix element of 
El matrix element obtained using DF method and RCC amplitude obtained using T^^^ operator. 
This shows that the upper most pi/2 core orbital, i.e. 6pif2 orbital, contributes almost through 
this term and the remaining core orbital contributions are small. This can be understood as the 
energy difference between the valence orbital, 7si/2 'ind 6pl/2 core orbital is small and hence it 
contributes large. 

As in the case of the core orbital contributions, we have also investigated the role of various 
virtual orbitals considering the most dominant operator, 5*2^, and El matrix element obtained 
using DF method. These results are reported in Table Hill As seen in this table, Tpi/2i 8Pi/2 ^'^'^ 
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TABLE III: Break down of DSj^ + c.c. contributions using DF reduced El matrix elements from important 
virtual orbitals. 

{lsi/2\\D\\npi/2)DF {npi/2\s''i^\7si/2) R 



n ~ 


7 


5.12 


227.26 


949.69 


n — 


8 


-0.57 


-81.71 


37.79 


n — 


9 


0.18 


115.55 


16.66 


n = 


10 


0.17 X 10-1 


-91.63 


-1.25 


n — 


11 


-0.53 X 10-1 


45.35 


-1.95 


n — 


12 


0.17 X 10-1 


-17.37 


-0.24 


n = 


13 


-0.21 X 10-2 


5.17 


-0.87 X 10 


n — 


14 


-0.29 X 10-3 


2.02 


-0.48 X 10 



9pi/2 orbitals contribute the most. In our DF calculations, the orbitals up to 8pi/2 are bound and 
9Pi/2 orbital onwards as continuum orbitals, respectively. Large contribution from the 7px/2 orbital 
can be understood based on the fact that the energy difference between the valence orbital 7si/2 
and virtual orbital 7pi/2 is very small and the overlap of these two orbitals s large at distances 
close to and far away from the nucleus. Since the density of the continuum orbital 9pi/2 in the 
nuclear region is large, it also gives a large contribution. 

We have also determined the contributions from the most important doubly excited states by 
combining the S^^l^ operator with the El matrix elements obtained using the DF method and given 
them in Table ITVl As a special condition, our 5*2^ can excite the valence orbital 7sx/2 to itself (a 
spectator), but other virtuals may go to the core orbitals by the dipole operator. This is a special 
type of core-polarization effect and which manifests through doubly excited states. From the above 
table, we find that virtual orbitals from the d symmetry are contribute significantly through this 
process. 

V. CONCLUSION 

We have discussed the relativistic theory of the electric dipole moments of paramagnetic atoms 
arising from the electric dipole moment of the electron. We have employed the relativistic coupled- 
cluster method to calculate the atomic wavefunctions to all orders in Coulomb interaction and 
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TABLE IV: Break down of DS^ J^ + c.c. contributions using DF reduced El matrix elements from important 
doubly excited states. Here indices m, I and k represent valence, core and virtual orbitals, respectively. 









{l\\D\\k)nF 


{ml\S^2\mk) 


R 


/ = 6pl/2; 


k = 


7s 1/2 


0.54 


11.05 


-5.25 


/ = 6pl/2; 


k = 


lOSi/2 


0.54 


13.78 


-2.22 


/ = 6pl/2; 


k = 


9^3/2 


-1.10 


-11.95 


-2.22 


/ = Qpl/2; 


k = 


10^3/2 


-1.48 


14.89 


-4.93 


I = 6p3/2; 


k = 


10^3/2 


-0.74 


-27.17 


-4.09 


/ = 6p3/2; 


k = 


7^5/2 


1.78 


-25.17 


-10.24 


/ = 6p3/2; 


k = 


8^5/2 


-1.18 


16.31 


-4.48 


I = 6p3/2; 


k = 


9^5/2 


-1.80 


24.87 


-10.55 


/ = 6p3/2; 


k = 


10^5/2 


-2.21 


27.31 


-14.80 


/ = 6p3/2; 


k = 


11^5/2 


1.10 


-7.85 


-2.50 



one order in the V and T violating interaction and obtained the EDM enhancement factor for 
francium; the heaviest of all the alkali atoms. Our result has been compared with the available 
semi-empirical result and contributions from various correlation effects and important intermediate 
states have been explicitly given. 
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